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Abstract: Vegetation near-soil-surface factors can protect topsoil from erosion, however, their 
contributions to the reduction of soil erosion, especially under natural rainfall events, have not been 
systematically recognized. This study was performed to quantify the effects of near-soil-surface factors on 
runoff and sediment under natural rainfall events on grasslands dominated by Bothriochloa ischaemum (Linn.) 
Keng (BI grassland) and Artemisia gmelinii Thunb. (AG grassland) in two typical watersheds on the Loess 
Plateau, China in 2018. By successive removal of the plant canopy, litter, biological soil crusts (BSCs) and 
plant roots, we established five treatments including plant roots, plant roots+BSCs, plant 
roots+BSCs+litter, intact grassland and bare land in each grassland type. In total, twenty runoff plots (5 
mX3 m) with similar slopes and aspects were constructed in the two types of grasslands. Results showed 
that plant canopy, litter and roots reduced runoff, while BSCs, which swelled in the presence of water, 
increased runoff. In contrast, all of these factors reduced sediment yield. In addition, the reductions in 
runoff and sediment yield increased with Iso (maximum 30-min rainfall intensity) for each vegetation 
near-soil-surface factor except for BSCs. Among these factors, plant canopy had the largest contribution 
to runoff reduction, accounting for 48.8% and 39.9% in the BI and AG grasslands, respectively. The 
contributions of these vegetation near-soil-surface factors to sediment yield reduction were similar 
(21.3%-29.9%) in the two types of grasslands except for BSCs in the AG grassland (10.3%). The total 
reduction in runoff in the BI grassland (70.8%) was greater than that in the AG grassland (53.1%), while 
the reduction in sediment yield was almost the same in both grasslands (97.4% and 96.7%). In conclusion, 
according to the effects of different vegetation near-soil-surface factors on runoff and sediment 
production, our results may provide more complete insight and scientific basis into the effects of various 
vegetation related factors in controlling soil erosion. 
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1 Introduction 


Soil erosion is a serious ecological problem worldwide. Many laboratory and field experiments 
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have confirmed that planting or preserving vegetation is an effective way to prevent soil loss, 
especially in vulnerable areas. Due to differences in vegetation types, vegetation distribution and 
community structure, the effects of vegetation on soil erosion vary (Lin et al., 2014; Chen et al., 
2018; Zhang et al., 2018). Compared with woodlands with a low coverage, grasslands with a high 
coverage or plant density may reduce annual runoff by 50% and annual soil loss by 92% (Chen et 
al., 2018). The various effects of vegetation on soil erosion are generally related to vegetation 
near-soil-surface factors such as plant canopy, litter, biological soil crusts (BSCs), root system and 
soil properties (Wang et al., 2014; Sun et al., 2016; Wang et al., 2017). 

Vegetation near-soil-surface factors include both aboveground and belowground factors, with 
the former including plant canopy, litter and BSCs, and the latter mainly including plant roots 
(Wang et al., 2014). Generally, plant canopy can intercept raindrops and reduce kinetic energy, 
consequently protecting the soil surface from splash erosion. Zhou et al. (2006) reported that for 
every 1% increase in plant canopy, the amount of soil loss decreased by 456 t/(km?.a) and soil 
erosion tended to be zero when the coverage reached 52%. Plant litter improves soil surface 
roughness and the combination of plant stems and litter forms a series of small dams, which 
decreases runoff rate and dissipates overland flow energy (Wang et al., 2014). Furthermore, plant 
litter can be incorporated into the topsoil and thus improve soil structure, enhance soil infiltration 
ability and promote soil stability (Liu et al., 2020). Plant litter incorporated into soil may decrease 
soil rill erodibility by 24.3%-34.8% in natural grasslands or woodlands (Sun et al., 2016). In 
contrast to plant canopy and litter, BSCs show both positive and negative effects on soil erosion. 
On the one hand, BSCs may promote soil water infiltration by increasing soil porosity and soil 
aggregate stability, thereby reducing overland flow velocity and erosive forces (Liu et al., 2016). 
On the other hand, BSCs may form an impermeable seal after absorbing water and thereby reduce 
infiltration and increase the overland flow rate (Coppola et al., 2011). Wang et al. (2014) reported 
that BSCs can be broken by overland flow and increase soil detachment capacity when the shear 
stress is greater than 11 Pa (i.e., the failure point). 

Plant roots are distributed throughout the soil mass and form a root-soil composite that 
strengthens soil shear strength during root growth (Gyssels et al., 2005). Plant roots also have 
strong direct effects on soil factors, increasing organic matter content, soil cohesion and 
water-stable aggregates and decreasing rill erodibility (Zhang et al., 2019). Wang et al. (2014) 
summarized the role of the root system in consolidating soil as involving two mechanisms: 
physical binding and chemical bonding. They showed that the root bonding effect accounted for a 
quarter of the total reducing effects of roots on soil erosion. In addition, root type greatly affects 
soil erosion. In general, plants with fibrous root systems are more effective in reducing soil 
erosion than those plants with tap roots (de Baets et al., 2007; de Baets et al., 2011). As reported 
by Wang and Zhang (2017), the ability of taproots to reduce soil detachment capacity is only 
two-thirds that of fibrous roots. In addition, soil properties (especially physical characteristics) 
have strong effects on soil erosion. An increase in soil bulk density reduces soil porosity and 
hinders infiltration, thus increasing overland flow and soil loss. Soil organic matter can enhance 
the formation of soil aggregates, improving soil structure and reducing soil erosion (Zhang et al., 
2019). 

However, most studies assessing the effects of vegetation on soil erosion have considered only 
one or two of the above factors. In addition, popular soil erosion prediction models ignore these 
parameters. For example, the cover-management factor has been used to reveal the effects of 
vegetation cover and management measures on soil erosion in the universal soil loss equation 
(USLE) and revised universal soil loss equation (RUSLE; Biesemans et al., 2000). However, the 
cover-management factor has been estimated at watershed and regional scales, mainly by the 
land-use valuation methods and vegetation coverage calculations using remote sensing images 
(Feng and Zhao, 2014). These methods do not consider the effects of some factors (such as plant 
roots or litter) on soil erosion, which may limit model precision. Moreover, it remains unknown 
the extent to which soil loss can be reduced by a single vegetation near-soil-surface factor and the 
contribution of the factor to the reduction of soil erosion under the participation of all of the 
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above-described vegetation near-soil-surface factors. Many studies have been carried out in this 
area and have achieved remarkable results. However, to our knowledge, only Wang et al. (2015) 
and Li and Pan (2018) took into consideration all vegetation near-soil-surface factors to determine 
the effects of vegetation on soil erosion. As reported by Wang et al. (2015), plant roots, BSCs and 
litter reduced the soil detachment capacity by 0.001—-0.094 kg/(m*-s), with contributions ranging 
from 17.1% to 58.7%. Li and Pan (2018) showed that plant roots were the primary contributor to 
reduced soil erosion, with an average contribution of 84%, while for runoff reduction, the 
contribution of aboveground vegetation near-soil-surface factors (plant canopy and litter) was 
almost 86%. However, they used either scouring modeling experiments or simulated rainfall on 
planted grassland (for two years) rather than natural rainfall on natural grassland. This may have 
caused their results to be unrepresentative of field conditions. A systematic study of all vegetation 
near-soil-surface factors (plant canopy, plant roots, BSCs and litter) under natural rainfall and 
natural grassland conditions may help us better understand their impacts on soil erosion and 
improve the prediction accuracy of soil erosion models. 

The Loess Plateau is among the areas with the most severe soil erosion in China. The Chinese 
government implemented a series of ecological construction projects to solve ecological problems 
in the Loess Plateau, such as the Grain for Green Project in 1999. Due to the implementation of 
this project, the vegetation coverage of the Loess Plateau increased from 31.6% in 1999 to 63.2% 
in 2018, and grassland became the main land-use type on the plateau, covering an area of 2.6x10° 
km? and accounting for 41.7% of the total plateau area in 2010 (Wang et al., 2015). The 
near-soil-surface factor of the grassland in the Loess Plateau is worth investigating due to the 20-a 
restoration. Wang et al. (2015) reported that plant litter accumulation in a grassland abandoned for 
24 a was 0.20 kg/m, which was 6.7 times that of grassland abandoned for 1 a. Plant litter reduced 
soil erosion by 8.1% and 13.3% in the grasslands abandoned for 1 and 24 a, respectively, 
compared with that of bare land (Wang et al., 2015). In addition, BSCs developed during 
vegetation restoration and covered up to 60%—70% of the soil surface on grassland of the hilly 
Loess Plateau. Furthermore, BSCs reduced sediment yield by 92% in revegetated grassland, and 
soil loss could be completely inhibited by BSCs in the later successional stage (Zhao and Xu, 
2013). Additionally, plant species composition and abundance have changed accordingly during 
natural succession (1—30 a) on grassland of the Loess Plateau (Sun et al., 2017). Different grass 
species with various root morphologies and structures, may have diverse effects on soil erosion. 
For example, in a previous study, the reduction in soil detachment capacity caused by Gramineae 
plants was 1.6-1.9 times that by Compositae plants due to the differences in their root systems 
(Wang and Zhang, 2017). Above all, previous studies concentrated on a limited near-soil-surface 
factors and failed to understand the overall effects of vegetation on soil erosion. 

In the present study, Bothriochloa ischaemum (Linn.) Keng and Artemisia gmelinii Thunb. are 
typical herbaceous plant species on grassland of the Loess Plateau, China. By observing natural 
rainfall events in the field, we studied the effects of vegetation near-soil-surface factors on soil 
erosion in these two types of grasslands. We hypothesized that all vegetation near-soil-surface 
factors reduce runoff and sediment, in which plant canopy performs better in reducing runoff 
whereas plant roots are most effective at decreasing sediment yield. The aims of this study were 
to: (1) investigate the reductions in runoff and sediment yield caused by vegetation 
near-soil-surface factors; (2) distinguish the relationships between runoff or sediment yield 
reductions and rainfall characteristics; and (3) quantify the relative contribution of each 
vegetation near-soil-surface factor to runoff and sediment yield reduction. 


2 Materials and methods 


2.1 Study area 


This study was carried out on two typical grasslands in two watersheds in Ansai District, Shaanxi 
Province in 2018. The grassland in the Fangta Watershed (8.4 km?; 36°47'27"N, 109°16'46”E) 
was dominated by A. gmelinii (abbreviated AG), and the grassland in the Xiannangou Watershed 
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(50.6 km?; 36°44'25"N, 109°17'22”E) was dominated by B. ischaemum (abbreviated BI; Fig. 1). 
The altitude of the study area ranges from 470 to 1800 m a.s.l. The loess-dominated hilly and 
gully region has a warm and temperate semi-arid climate. The annual mean temperature is 8.8°C, 
and the mean annual precipitation is 400-500 mm. Heavy rainfall events often occur in summer 
(July-September) and result in severe soil erosion. This region belongs to the forest grassland 
transition zone, and Compositae and Gramineae are typical herbaceous plant taxa in this area. The 
original land-use type of both two types of grassland was sloping farmland, with a slope gradient 
of 60%, and both grassland areas have experienced natural succession for approximately 35 a 
since they were abandoned. The physical properties of the test soil are shown in Table S1. 
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Fig. 1 (a), distribution of two experimental plots; (b), BI grassland is the grassland with B. ischaemum as the 
dominant species; (c), AG grassland is the grassland with A. gmelinii as the dominant species. A schematic 
diagram of T1—T4 treatments is shown in Figure S1. The treatments named TO, T1, T2, T3, T4 and T5 indicate 
bare land, plant roots, plant roots+BSCs, plant roots+BSCs+litter, intact grassland and dead roots, respectively, 
and the T5 is not analyzed in the study. BSCs, biological soil crusts. 
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2.2 Experimental design 


Four treatments including bare land (TO), plant roots (T1), plant rootstBSCs (T2), plant 
roots+BSCs+litter (T3) and intact grassland (T4) were set up to investigate and distinguish 
different vegetation near-soil-surface factors in controlling soil erosion. For bare land (TO), the 
0-30 cm topsoil was flipped, sieved and packed into the plots to a bulk density, which was used 
as the control. Two runoff plots were designed for each treatment as replicates, and a total of ten 
runoff plots were constructed on each type of grassland in June 2018. The runoff plots were 
surrounded by 40-cm-high PVC sheets buried in the soil to a depth of 20 cm. Each runoff plot had 
an area of 15 m? (5 m in length and 3 m in width). 


2.3 Rainfall characteristics 


The rainfall amount (P) and duration (t) were recorded by a tipping bucket self-recording rainfall 
gauge (RG3-M, HOBO, USA). During the observation period (July-September in 2018), a total 
of 20 and 15 rainfall events were detected by monitoring equipment in the BI and AG grasslands, 
respectively (Fig. S2). Among them, 10 and 7 rainfall events produced runoff and sediment, 
respectively, and these events were termed effective rainfall events. To eliminate the disturbance 
caused by the building of runoff plots, we did not collect the first runoff and sediment samples. 
Therefore, a total of 15 effective rainfall events (9 rainfall events in the BI grassland and 6 in the 
AG grassland) were analyzed in this study (Table S2). After each rainfall event, P and t data were 
collected and used as the basic data for the next analysis. Other rainfall parameters, namely, 
average rainfall intensity (Za), maximum 30-min rainfall intensity (/30), rainfall kinetic energy (E) 
and rainfall erosivity (R) were calculated to describe the rainfall characteristics using Rain Record 
software. 


2.4 Measurements of runoff and sediment yield 


For each plot, the runoff and sediment yield were collected after each rainfall event by using a 
tank that was fixed to the bottom of runoff plot. The volume of runoff-sediment mixture in the 
tank (Vsr-p) was measured with a ruler. Next, the muddy water in each tank was fully mixed to 
ensure that it was uniform and sampled with a 1000-mL plastic bottle, and this procedure was 
repeated twice. The volume of muddy water in each plastic bottle was measured (Vsr-») and 
weighed (Msr-»). Then, the muddy water was moved into an aluminum box, the sediment was 
allowed to settle, the supernatant was drained and the sediment was dried (105°C, 24 h) to obtain 
the weight of sediment yield (Ms-p) in the plastic bottle. The ratio of sediment to runoff was used 
to calculate the amount of runoff (Vr-p) and the sediment yield (Ms-p) of each runoff plot. 


M,, = Msp» Vrp Ms, ; (1) 
Vern Msp» 
Vag =e ee, (2) 


where Msr» and Vsr-» are the weight (g) and volume (mL), respectively, of muddy water in the 
plastic bottles; Msr-p and Vsr-p are the weight (g) and volume (mL), respectively, of muddy water 
in the runoff plot; Ms» and Ms-p are the weight (g) of sediment yield in the plastic bottle and 
runoff plot, respectively; and pw is the density of clean water (g/cm?). 

In this study, runoff generation was expressed by runoff depth (RD, mm), and soil loss was 
measured as sediment yield per unit area (SY, g/m7), calculated by the following equations: 


RD= VR, (3) 
A 

sy =- ss, (4) 
A 


where A is the projected area of the runoff plot (15 m7). 
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2.5 Contributions of near-soil-surface factors to soil erosion reduction 


In this study, the contribution to runoff or sediment reduction was applied to represent the 
inhibitory effects of near-soil-surface factors on soil erosion. The calculation process was as 
follows: The participating near-soil-surface factors differed among the treatments, and the runoff 
and sediment yield varied. Therefore, the reductions in runoff and sediment yield caused by 
near-soil-surface factors of the plant canopy, plant roots, BSCs and litter were calculated by 
Equations 5 and 6. 
AR; = Ry — Ry (5) 
AS p = Sp Sra (6) 
where R is the runoff amount of each treatment (mm); S is the sediment yield (g/m7); and T; is the 
treatment, with i=1, 2, 3, 4. 
The actual contribution of each factor i (plant canopy, plant roots, BSCs and litter) to reducing 
runoff (CR, %) or sediment (CS, %) was computed by Equations 7 and 8. 


cR,, = Ka x100% (7) 


Ti 
CS, = Sr 100% . (8) 
Sri 
The actual contributions of T4 (intact grassland) to runoff (CR74, %) and sediment yield 
(CSz, %) reductions were expressed by Equations 9 and 10. 


CRra — Bro — rs 100%, (9) 
Rra 

CS _ Sro Sra 199%. (10) 
S 


T4 
2.6 Statistical analysis 


Because all runoff plots on the same grassland were subjected to the same rainfall events, the 
reduction in runoff or sediment yield or the contribution to runoff or sediment yield reduction 
could be compared among different treatments within a grassland. We calculated runoff and 
sediment reduction amounts and the corresponding contributions of each near-soil-surface factor 
under each effective rainfall event according to the above formulae, and used the results as the 
basic data for the following analysis. Spearman correlation analysis was used to evaluate the 
relationships between rainfall parameters and the contributions of near-soil-surface factors to 
runoff and sediment reduction. A curve estimation model was used to identify the relationship 
between the amount of runoff or sediment reduction and 730. Additionally, Kruskal-Wallis analysis 
was used to test for differences in the contribution to runoff or sediment yield reduction among all 
of the near-soil-surface factors. All analyses were performed using SPSS v19.0 software, and 
graphs were constructed with Origin 2018 software. 


3 Results 
3.1 Runoff and sediment yield 


The runoff and sediment yield varied greatly among treatments, with near-soil-surface factors 
successively superimposed (corresponding to T1—T4) (Table S3 and S4). Specifically, the runoff 
yield ranged from 0.0 to 5.3 mm in the BI grassland and from 0.0 to 7.0 mm in the AG grassland. 
Compared with that in bare land, runoff was reduced by 39.7%, 43.9%, 75.1% and 93.4% in T1, 
T2, T3 and T4 treatments, respectively, in the BI grassland, and the corresponding values in the 
AG grassland were 11.4%, 7.3%, 47.4% and 72.2%, respectively. The ranges of sediment yield in 
the BI and AG grasslands were 0.0-1138.8 and 0.0-476.8 g/m’, respectively, with the average 
values of 53.0 and 32.1 g/m?, respectively. The different treatments exerted different effects on 
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sediment yield. For example, T1, T2, T3 and T4 reduced sediment yield relative to that of bare 
land by 65.1%, 92.4%, 98.7% and 99.7%, respectively, in the BI grassland and by 57.2%, 76.4%, 
93.6% and 96.9%, respectively, in the AG grassland. 


3.2 Reduction of vegetation near-soil-surface factors on runoff and sediment yield 


Overall, the effects on runoff and sediment reduction differed among the vegetation near-soil- 
surface factors, i.e., plant canopy, litter, plant roots and BSCs, in the two types of grasslands (Fig. 
2). The plant canopy, litter and plant roots contributed to runoff reduction in both BI and AG 
grasslands, with litter being the greatest contributor (Fig. 2a and 2b). The amount of runoff 
reduction achieved by plant canopy, litter and plant roots was 0.3, 0.6 and 0.5 mm, respectively, in 
the BI grassland and 0.4, 0.7 and 0.2 mm, respectively, in the AG grassland. BSCs showed an 
opposite effect compared with other factors, contributing to runoff yield. 

Similarly, all near-soil-surface factors reduced sediment yield. In the BI grassland, plant canopy, 
litter, BSCs and plant roots reduced the sediment yield by 0.0-649.9 g/m? (Fig. 2d). Among the 
factors, plant roots achieved the highest average sediment reduction (108.0 g/m), which was 56.0, 
7.9 and 2.1 times that of plant canopy, litter and BSCs, respectively. In the AG grassland, the 
near-soil-surface factors reduced sediment yield by 0.0-282.5 g/m? (Fig. 2d). The plant roots 
produced the greatest reduction in sediment yield (53.8 g/m?), which were 17.0, 3.3 and 2.9 times 
those of plant canopy, litter and BSCs, respectively. 
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Fig.2 Reduction of runoff amount and sediment yield of different near-soil-surface factors at each rainfall event 
in the BI grassland (a, c) and AG grassland (b, d). F1, F2, F3 and F4 refer to plant canopy, litter, BSCs and plant 
roots, respectively. 
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3.3 Contributions of vegetation near-soil-surface factors to reductions of runoff and 
sediment yield 


In general, the total contribution of all of the near-soil-surface factors to reducing runoff in the BI 
grassland was 70.8% (ranging from 19.5% to 94.7%), which was 1.3 times greater than that in the 
AG grassland (53.1%; ranging from —19.4% to 79.4%). In addition, the factors in both BI and AG 
grasslands made significant contributions to controlling sediment loss, which was reduced by 
97.4% (ranging from 92.2% to 99.8%) in the BI grassland and 96.7% (ranging from 92.9% to 
99.1%) in the AG grassland. 

In the BI grassland, the contribution of plant canopy to runoff reduction was 48.8%, which was 
3.6 and 3.3 times those of litter (13.5%) and plant roots (14.8%), respectively (Table 1). In the AG 
grassland, the contribution of plant canopy was 39.9%, which were 1.7 and 7.5 times those of 
litter (23.9%) and plant roots (5.3%) (Table 1). For BSCs, their contributions to reducing runoff 
were negative, being —6.4% (algae) in the BI grassland and —16.1% (moss) in the AG grassland, 
respectively. The contributions of factors (except BSCs) to soil sediment reduction were similar, 
with percentage of approximately 25% (Table 1). Plant canopy, litter, BSCs and plant roots made 
contributions of 21.3%, 24.3%, 24.2% and 27.7%, respectively, in the BI grassland and 29.5%, 
29.9%, 10.3% and 27.0%, respectively, in the AG grassland. However, the contribution of BSCs 
to sediment yield reduction in the AG grassland was only 10.3%, which was lower than the 
contributions of other factors. 


Table 1 Average contribution of each near-soil-surface factor to runoff and sediment yield reduction 


Near-soil-surface Runoff Sediment yield 

factor BI grassland AG grassland BI grassland AG grassland 
Plant canopy (%) 48.8° 39.9% 21.3? 29.5° 
Plant litter (%) 13.5° 23.9" 24.3* 29.9° 
BSCs (%) —6.4° —16.1° 24.27 10.3° 
Plant roots (%) 14.8° 5.3° 27.7" 27.08 


Note: Different lowercase letters within a column indicate significant differences at P<0.05 level. BI grassland is the grassland with B. 
ischaemum as the dominant species, and AG grassland is the grassland with A. gmelinii as the dominant species. BSCs, biological soil 
crusts. 


4 Discussion 


4.1 Effects of vegetation near-soil-surface factors on reduction of runoff and sediment 
yield 

With the gradual improvement of vertical structure of vegetation community in the grassland, the 
reductions in runoff and sediment yield gradually strengthened. Wen et al. (2010) found that a 
complex stratified structure was more effective than a single vegetation layer in controlling soil 
erosion. For example, in contrast to single factors, multiple near-soil-surface factors reduced 
runoff generation by 51.0% and sediment yield by 72.3% due to their strong ability to absorb 
raindrop energy through multiple interceptions (for example, plant canopy and litter) (Zheng et al., 
2008). Our results showed that runoff and sediment in both BI and AG grasslands were reduced 
by 7.3% to 93.4% compared with that in bare land (Tables S3 and S4), which represented simple 
near-soil-surface factors. Furthermore, we confirmed that different near-soil-surface factors 
exerted different effects on runoff and sediment yield. 

In the BI and AG grasslands, plant litter was more effective in reducing runoff and sediment 
yield than plant canopy, plant roots and BSCs. This result can be attributed to various water 
holding capacity abilities of different factors. Generally, plant litter reduces overland flow by 
covering the ground and absorbing water (Facelli and Pickett, 1991). For example, water holding 
capacity of plant litter in typical plant communities was found in the range from 5.1 to 12.8 g/g in 
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the hilly region of the Loess Plateau (Wang et al., 2018). In addition, plant litter can be 
incorporated into the topsoil and hold water, thereby reducing runoff production (Sun et al., 2016). 
More importantly, plant litter covering the soil surface improves soil surface roughness and 
reduces flow velocity and energy of overland flow, thereby inhibiting soil loss (Wang et al., 2015). 
Furthermore, plant litter incorporated into the topsoil can improve soil structure, permeability and 
stability, thus reducing soil detachment capacity (Sun et al., 2016). In a previous study, the 
incorporation of grass litter into topsoil (with a litter rate of 0.2 kg/m? and a soil bulk density of 
1.25 g/cm?, respectively) reduced soil loss compared with that on bare land by 34% and 21% 
under rainfall intensities of 40 and 80 mm/h, respectively (Wang et al., 2020). Moreover, plant 
canopy was found to reduce runoff production by intercepting raindrops and delay the time of 
runoff initiation (Zhou et al., 2006), and plant roots systems might improve soil structure (e.g., 
soil porosity) and enhance soil water infiltration capability (Gyssels et al., 2005). 

However, our study showed that BSCs increased runoff, by an average of 0.0 mm (ranging 
from —0.5 to 0.2 mm) in the BI grassland and an average of 0.1 mm (ranging from —19.4 to 0.3 
mm) in the AG grassland (Fig. 2a and 2b). Although several studies have reported that the growth 
of BSCs promotes the infiltration capacity of soil by increasing soil porosity and enhancing soil 
aggregate stability, thereby reducing runoff (Zhao and Xu, 2013). BSCs can form an impermeable 
seal by absorbing water during rainfall events, which may inhibit soil water infiltration and 
accelerate runoff generation (Coppola et al., 2011). In the present study, BSCs increased runoff 
but reduced sediment yield. The reasons may be as follows: (1) BSCs increased soil surface 
roughness, reduced flow velocity and weakened flow shear stress (Kheirfam and Roohi, 2020); (2) 
rainfall intensity was low, but rainfall duration was long, resulting in low flow velocity and 
energy; and (3) flow shear stress was not enough to destroy BSCs, and consequently, sediment 
yield did not increase with increasing runoff (Wang et al., 2014). Overall, our results indicate that 
near-soil-surface factors play essential roles in reducing runoff and sediment and that different 
factors show diverse effects. 


4.2 Contributions of vegetation near-soil-surface factors to runoff and sediment yield 
reduction 


Our results showed that the total contribution of all of the near-soil-surface factors to reducing 
runoff in the BI grassland was 70.8% (ranging from19.5% to 94.7%), which was 1.3 times greater 
than the corresponding contribution in the AG grassland (53.1%; ranging from —19.4% to 79.4%). 
This difference can be attributed to the different structural and morphological traits of the two 
dominant herbaceous plants in the two grasslands. In the BI grassland, the canopy has a 
tussock-forming architecture and covers the soil surface like a carpet, which is effective in 
intercepting rainfall. In the AG grassland, plant canopy features an expanding shoot architecture, 
and plants are spaced far apart, reducing the chance of intercepting raindrops (Du et al., 2013). 
Plant litter in the BI grassland is standing litter and is not effective in reducing runoff. However, 
plant litter in the AG grassland is composed of thick branches and small leaves scattered on the 
soil surface, which can hinder and accumulate runoff. The fibrous root systems in the BI 
grassland have numerous fine roots, which increase soil porosity and water infiltration more 
effectively than the tap root systems in the AG grassland. All these differences resulted in lower 
runoff reduction in the AG grassland than in the BI grassland. 

Because near-soil-surface factors affected runoff in different ways and via different 
mechanisms, the contributions of plant canopy, litter, BSCs and plant roots to runoff reduction or 
increase differed. Among all of the near-soil-surface factors, plant canopy had the maximum 
contributions to runoff reduction in both BI and AG grasslands. In the BI grassland, the 
contribution of plant canopy to runoff reduction was 48.8%, which was 3.6 and 3.3 times those of 
litter (13.5%) and plant roots (14.8%), respectively (Table 1). In the AG grassland, the 
contribution of plant canopy was 39.9%, which was 1.7 and 7.5 times those of litter (23.9%) and 
plant roots (5.3%; Table 1). These differing contributions to runoff reduction were probably 


JOURNAL OF ARID LAND 


related to different spatial distribution of these factors. During rainfall events, plant canopy is the 
first structure to intercept rainfall, and its effects on runoff are not affected by other factors (e.g., 
litter and plant roots). Regarding BSCs, their contributions to runoff reduction were negative, i.e., 
—6.4% (algae) in the BI grassland and —16.1% (moss) in the AG grassland, respectively (Table 1). 
BSCs increased runoff generation through the formation of an impermeable layer after swelling, 
and the swelling mechanisms of different components of BSCs are different. As reported by Liu et 
al. (2016), through the hydration of the sheath, the expansion of algae occurs, while that of moss 
occurs through the hydration of stems and leaves. Furthermore, the swelling height of moss was 6 
times that of algae, which resulted in moss being more impermeable and increasing runoff 
generation more than algae. 

The BI and AG grasslands both made significant contributions to the control of soil erosion, 
which was reduced by 97.4% (ranging from 92.2% to 99.8%) in the BI grassland and 96.7% 
(ranging from 92.9% to 99.1%) in the AG grassland (Table 1). These results indicate that 
long-term restored grasslands have a strong ability to control soil erosion. Differences in the 
reductions in sediment yield among near-soil-surface factors were not significant despite the 
variation in the mechanisms controlling soil erosion, as mentioned above. In other words, the 
contributions of these factors (except for BSCs) to sediment reduction were similar, with 
proportions of approximately 25% (Table 1). Plant canopy, litter, BSCs and plant roots made 
contributions of 21.3%, 24.3%, 24.2% and 27.7%, respectively, in the BI grassland and 29.5%, 
29.9%, 10.3% and 27.0%, respectively, in the AG grassland (Table 1). Therefore, 
near-soil-surface factors all contributed to the control of soil loss. However, the contribution of 
BSCs to sediment yield reduction in the AG grassland was only 10.3%, which was lower than the 
contributions of other factors. As mentioned above, BSCs can reduce soil infiltration and enhance 
runoff generation; thus, they might be expected to exacerbate soil loss. In the AG grassland, BSCs 
were composed mainly of moss, which have a strong ability to promote runoff production. Wang 
et al. (2014) reported that BSCs accelerated soil erosion after being fractured when the flow shear 
stress exceeded its threshold. 


4.3 Effects of rainfall characteristics on soil erosion mediated by near-soil-surface factors 


Rainfall is an essential factor affecting runoff and sediment yield on slopes, and it can be 
described by various parameters, such as P, Ia, t, Bo, E and R. In general, P, t and Ia are 
commonly used to reflect the effects of rainfall on soil erosion, and generally, runoff and sediment 
yield increase with increases in these parameters. In this study, among all 6 parameters 
representing rainfall characteristics, there was a strong correlation between /30 and each of other 
rainfall parameters as well as the reduction in runoff or sediment yield (Fig. 3). In addition, as 
mentioned by Renard (1997), Iso can represent the difference in rainfall intensity and R, and 
reflect the prolonged peak rates of detachment and runoff during a whole rainfall event. Therefore, 
Bso was applied in this study to reflect the effects of rainfall characteristics on soil erosion under 
all rainfall events. 

In this study, the reduction in runoff or sediment yield of each near-soil-surface factor (with the 
exception of BSCs in the BI grassland) increased with Jz0 (Fig. 4). In the BI grassland, the 
reductions in runoff and sediment yield caused by plant canopy, litter and plant roots increased 
with Jo as a power function (R7>0.74; P<0.01; Fig. 4). In the AG grassland, only runoff 
reductions by plant canopy and litter increased significantly with J30 as a power function (R7>0.75; 
P<0.01), and the reduction in runoff caused by plant roots and the reduction in sediment yield 
caused by plant canopy, litter and plant roots tended to increase with increasing Zso (Fig. 4). In 
addition, when /39 was greater than 10 mm/h, the runoff and sediment reductions of plant canopy, 
litter and plant roots in both BI and AG grasslands increased rapidly with increasing /39. Moreover, 
when taking intact grassland into consideration, power functions were detected between rainfall 
characteristics (represented by J30) and total contributions to runoff and sediment yield reductions 
in the BI and AG grasslands (P<0.05; Fig. 5). As Io increased, the reduction in runoff or sediment 
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yield increased continuously (Fig. 5). This pattern could have been observed because the ability 
of all of the near-soil-surface factors to reduce runoff or sediment yield was not fully reflected in 
natural rainfall events (with /30 being less than 12 mm/h in this study); moreover, vegetation 
influenced the reductions in runoff and sediment yield to a greater extent than rainfall conditions. 
However, when /39 reaches a certain value, the ability of near-soil-surface factors to inhibit soil 
erosion might reach a limit. Therefore, it is necessary to further explore the rainfall conditions 
under which contribution of near-soil-surface factors to the inhibition of soil erosion are the 
greatest and stable. 
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Fig. 3 Correlation analysis of rainfall characteristics with runoff and sediment yield under different 
near-soil-surface factors. P, precipitation; t, rainfall duration; Ja, average rainfall intensity; /30, maximum 30-min 
rainfall intensity; E, rainfall kinetic energy; R, rainfall erosivity. RR_C and SR_C indicate runoff and sediment 
yield reduction caused by plant coverage; RR_L and SR_L indicate runoff and sediment yield reduction caused by 
litter; RR_B and SR_B indicate runoff and sediment yield reduction caused by BSCs (biological soil crusts); and 
RR_R and SR_R indicate runoff and sediment yield reduction caused by plant roots. **, P<0.01 level;*, P<0.05 
level. 
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Fig. 4 Variations in runoff and sediment yield reduction with /30 (maximum 30-min rainfall intensity) under 
different near-soil-surface factors in the BI grassland and the AG grassland. (a, e), plant canopy; (b, f), litter; (c, g), 
BSCs (biological soil crusts); (d, h), plant roots. 
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Fig. 5 Contributions of grassland to runoff reduction (a) and sediment yield reduction (b) 


5 Conclusions 


This study was performed to investigate the effects of near-soil-surface factors (plant canopy, 
litter, BSCs and plant roots) on runoff and sediment yield under natural rainfall events by 
successive removal of these factors in field. Our results showed that plant canopy, litter and plant 
roots reduced the runoff amount, while BSCs promoted runoff generation. All these 
near-soil-surface factors reduced the sediment yield. The reductions in runoff and sediment yield 
caused by plant canopy, litter and plant roots either tended to increase with 73o or increased with 
Bo as a power function. Among all of the near-soil-surface factors, plant canopy had the largest 
contribution to runoff reduction, which was much greater than those of other factors, in both two 
types of grasslands, and BSCs contributed the most to runoff production. All of the 
near-soil-surface factors contributed equally to reducing sediment yield (by approximately 25%) 
except for the moss crusts in the grassland dominated by A. gmelinii (approximately 10.3%). In 
general, grass is more effective at reducing soil loss than runoff. Moreover, due to the differences 
in structural composition and morphological characteristics between two grassland communities, 
the grassland with B. ischaemum as the dominant species performed better at reducing runoff. 

Due to the unique climate and vegetation conditions in arid and semi-arid regions, it is essential 
and practical to promote natural restoration of grasslands in these regions to prevent soil erosion 
in the long term. This study enhances our understanding of the mechanisms by which grassland 
communities conserve soil and water, which can improve the accuracy of many soil erosion 
models, especially those for arid and semi-arid regions. This study also has reference value for 
assessing actual soil and water loss situations and guiding the implementation of ecological 
projects. 
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Appendix 
Table S1 Basic soil information of sample plots 
Sample plot Soil layer Soil organic Total nitrogen NOs-N NH4*-N Bulk density 
site (cm) carbon (g/kg) (g/kg) (mg/kg) (mg/kg) (g/cm?) 
0-10 6.33 0.39 3.16 3.98 1.22 
BI 10-20 3.09 0.21 0.79 2.71 1.41 
grassland 20-30 2.68 0.18 0.77 4.79 1.46 
30—40 2.50 0.16 0.58 2.59 1.44 
0-10 7.67 0.41 0.55 3.64 1.16 
AG 10-20 4.19 0.29 0.82 2.67 1.19 
grassland 20-30 3.17 0.18 0.66 3.18 1.33 
30—40 2.76 0.17 0.63 3.44 1.31 


Sample plot Soil layer Soil water-stable aggregate (%) 


site (cm) >5.00 mm 2.00-5.00 mm 1.00-2.00 mm 0.50-1.00 mm 0.25-0.50 mm <0.25 mm 
Meot 0-20 33.96 8.47 6.83 6.60 5.05 39.09 
AG 
grassland 0-20 21.35 6.94 4.19 4.00 2.94 60.58 


Note: BI grassland is the grassland with B. ischaemum as the dominant species, and AG grassland is the grassland with A. gmelinii as the 
dominant species. 


Table S2 Characteristics of rainfall events 


Sample site Rainfall date Rainfall Rainfall duration Bo Rainfall kinetic Rainfall GrosiVily 
(mm) (min) (mm/h) energy (MJ/hm^) (MJ-mm/(hm*-h)) 
31 Jul 5.6 684 7.3 1.0 7.3 
9 Aug 3.6 263 4.8 0.5 2.5 
11 Aug 3.0 113 5.6 0.6 3.2 
18 Aug 9.8 147 12.0 2.0 23.9 
Pe 22 Aug 41.0 1383 11.6 63 73.1 
31 Aug 15.8 890 5.8 2.2 12.7 
2 Sep 16.0 351 10.6 ZT 28.5 
14 Sep 15.8 542 11.6 2.5 28.6 
19 Sep 13.2 3181 2.5 1.4 3.4 
9 Aug 4.4 227 6.8 0.9 6.3 
11 Aug 7.8 122 6.8 1.7 11.2 
AG 22 Aug 15.6 363 9.6 2.9 27.5 
grassland 31 Aug 26.4 779 11.6 5.1 59.5 
2 Sep 20.4 283 11.6 4.1 47.5 
14 Sep 9.6 498 5.0 1.7 8.2 


Notes: /30, maximum 30-min rainfall intensity. BI grassland is the grassland with B. ischaemum as the dominant species, and AG 
grassland is the grassland with A. gmelinii as the dominant species. 
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Table S3 Variations in runoff and sediment yield in the BI grassland under rainfall events 


Treatment Date pe pe l Date P ie Date a j 
TO 0.1 6.1 5.3 1138.8 1.4 38.8 
Tl 0.1 1.8 3:2 416.0 1.0 13.4 
T2 31 Jul 0.1 0.6 18 Aug 2.8 89.8 2 Sep 0.7 2.9 
T3 0.1 0.6 1.3 15.4 0.3 0.5 
T4 0.0 0.1 0.3 3.9 0.1 0.2 
TO 0.0 0.4 5.1 359.4 2.7 102.8 
Tl 0.0 0.2 3.3 118.3 1.4 25.8 
T2 9 Aug 0.0 0.1 22 Aug 3:3 22.2 14 Sep 1.3 8.8 
T3 0.0 0.0 1.3 3.2 0.6 1.2 
T4 0.0 0.0 0.3 0.6 0.1 0.3 
TO 0.0 1.1 0.3 9.1 0.0 0.1 
T1 0.0 0.3 0.2 2.4 0.0 0.0 
T2 11 Aug 0.0 0.2 31 Aug 0.2 0.6 19 Sep 0.0 0.0 
T3 0.0 0.1 0.1 0.2 0.0 0.0 
T4 0.0 0.0 0.1 0.2 0.0 0.0 


Notes: TO, bare land; T1, plant roots; T2, plant roots+BSCs; T3, plant roots+BSCs+ litter; T4, intact grassland. BSCs, biological soil 
crusts. BI grassland refers to the grassland with B. ischaemum as the dominant species. 


Table S4 Variations in runoff and sediment yield in the AG grassland under rainfall events 


Treatment Date pues ae t Date a. a 
TO 0.0 0.6 7.0 476.8 
Tl 0.0 0.1 3:5 202.3 
T2 9 Aug 0.0 0.1 31 Aug 5.7 104.8 
T3 0.0 0.0 3.3 30.1 
T4 0.0 0.0 2.1 15.8 
TO 0.0 0.5 2.3 46.0 
Tl 0.0 0.3 2.9 21.8 
T2 11 Aug 0.0 0.3 2 Sep 2.7 19.2 
T3 0.0 0.1 1.6 4.6 
T4 0.0 0.0 0.5 1.0 
TO 0.6 23.9 0.0 0.1 
Tl 0.4 9.9 0.0 0.1 
T2 22 Aug 0.7 4.9 14 Sep 0.0 0.1 
T3 0.3 0.3 0.0 0.1 
T4 0.1 0.2 0.0 0.0 


Note: TO, bare land; T1, plant roots; T2, plant roots+BSCs; T3, plant roots+BSCs+ litter; T4, intact grassland. BSCs, biological 
soil crusts. AG grassland refers to the grassland with A. gmelinii as the dominant species. 
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T4: Intact grassland 


T3: Plant roots 
+BSCst+litter 


T2: Plant roots 
+BSCs 


T1: Plant roots 


Fig. S1 Schematic diagram of vegetation near-soil-surface factors involved in different treatments (T1—-T4). 
BSCs, biological soil crusts. 
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Fig. S2 Recorded rainfall events in the BI and AG grasslands during rainy season. BI grassland is the grassland 
with B. ischaemum as the dominant species, and AG grassland is the grassland with A. gmelinii as the dominant 
species. Non-effective rainfall events indicate that no runoff or sediment is generated. Effective rainfall events 


indicate that runoff and sediment are generated. 


